Introduction

A conventional transistor is a three terminal solid-state semiconductor device which can be used for several applications in electronics. The basic use of a transistor is as a switch but with the technology shrinking, the researchers have started working to produce the smaller and faster switch. As stated by Moore’s law, for every 2 years the number of transistors on a chip should be doubled but as the transistor size is entering into the nanometer regime, transistors are suffering from several second-order effects and dominance of leakage current (Mollick 2006). To meet out these challenges, the single-electron transistors (SET) are likely to become an alternative for producing smaller and faster switches. The ultra small size of SET will allow the designers to integrate more number of transistors on a chip. In recent years, researchers have started focusing on the use of the organic molecules as island in the single-electron transistor (Stokbro 2010; Brandbyge et al. 2002; Neaton et al. 2006; Seldenthuis et al. 2010; Parashar et al. 2012).

Single-electron transistor (SET): Theory

Single-electron transistor is a device based on the concept of one-by-one electron tunneling through island utilizing the coulomb blockade effect. It consists of two electrodes named as drain and source, connected through tunnel junctions to a common electrode with a low self-capacitance, known as the island. The electrical potential of the island can be tuned by a third electrode, known as the gate which is capacitively coupled to the island as shown in Fig. 1. The only way to an electron to travel from one electrode to another is tunneling through the island. (Kumar and Kaur 2010; Fulton and Dolan 1987; Stokbro 2010; Grabert and Devoret 1992).

Fig. 1
figure 1

Schematic of single electron transistor

For an electron to tunnel onto the island, its energy must be equal to the Coulomb energy e2/2C, where C is the self-capacitance of the island defined as

$$ C = \frac{{e^{2} }}{\Updelta E} $$

Here, \( \Updelta E \) is the separation between the energy levels of the island.

Capacitance of the island must be <10−17 Farads and, therefore, its size should not be >10 nm. In SET at zero gate and bias voltages, the electrons do not have enough energy to enter the island and hence, there will be no current flow. When the bias voltage between the source and drain increases, electron can pass through the island as its energy reaches the Coulomb energy. This effect is known as the Coulomb blockade. The critical voltage needed to transfer an electron onto the island, equal to e/C, is called as Coulomb gap voltage. To achieve the Coulomb blockade, three criteria have to be met:

  1. 1.

    The bias voltage must be lower than the elementary charge divided by the self-capacitance of the island: \( V_{bias} < \frac{e}{C} .\)

  2. 2.

    The thermal energy in the source contact plus the thermal energy in the island i.e., K B T, must be below the charging energy: \( K_{B} T < \frac{{e^{2} }}{C}, \) or else the electron will be able to pass the island via thermal excitation.

  3. 3.

    And the tunneling resistance, Rt, should be greater than \( \frac{h}{{e^{2} }} \) which is derived from Heisenberg’s uncertainty principle.

The Coulomb blockade tries to alleviate any leak by current during the off state of the single-electron transistor (Fulton and Dolan 1987; Grabert and Devoret 1992; Altshuler et al. 1991).

In SET environment, electron flow happens because of the difference of agenda of Fermi functions in the two contacts where one keeps filling up the level(s) while the other empties them, the result of this process is a net flow of current. Here, all the levels do not conduct, only some levels support the flow of electrons. From Fig. 2, one can see that the two Fermi functions that levels deep down and way above do not contribute to current, because at those levels f1 and f2 have the same values of 1 and 0’s. The maximum conductance of SET is also called as the Quantum of conductance, G max  = 2q2/h. Correspondingly, the minimum resistance value is R min  = h/2q2.

Fig. 2
figure 2

Electron tunneling through the levels

The current in the single-electron transistor and the average number of electrons can be given by the expressions:

$$ I = \frac{q}{h} \frac{{\gamma_{1} \gamma_{2} }}{{\gamma_{1 + } \gamma_{2} }} (f_{1} - f_{2} ) $$
(1)
$$ N = \frac{{\gamma_{1} f_{1} + \gamma_{2} f_{2} }}{{\gamma_{1} + \gamma_{2} }} $$
(2)

Here \( \gamma_{1 } \) and \( \gamma_{2 } \) are the parameters which explains the electron escaping from source or island. On dividing these parameters by \( h \), the rate of escape can be explained. The last factor in the current expression tells that current flows if there is a difference between the Fermi functions of contacts. If both f1 and f2 are equal, then current will be zero (Datta 2004).

The validity of above expressions lies with the fact that applied voltage must be small. If the coupling is stronger, then we have to take into account the concept of ‘Broadening’. If a strong voltage is applied to the contacts, then the sharp discrete level broadens out and becomes a continuum of states as showed in Fig. 3, the sum of which is still one level and can accommodate one electron. In Fig. 3, if the energy range dE has density of states D(E), then the expressions for current and the average number of electrons given in the Eq.(1) and Eq.(2) can be expressed as, (Datta 2004).

$$ I = \frac{ - q}{h} \mathop \int \nolimits D\left( E \right)dE \frac{{\gamma_{1} \gamma_{2} }}{{\gamma_{1 } + \gamma_{2} }}(f_{1} - f_{2} ) $$
(3)
$$ N = \mathop \int \nolimits D\left( E \right)dE \frac{{\gamma_{1} f_{1} + \gamma_{2} f_{2} }}{{\gamma_{1} + \gamma_{2} }} $$
(4)
Fig. 3
figure 3

Broadening of island levels

To estimate the potential energy of the island, SET is represented through the model shown in Fig. 4, where CS, CD, and CG represent the capacitances established by the source, drain tunnel junctions, and gate dielectric. The amount of charge in the island is given by,

$$ Q = - q\Updelta n = C_{\text{S}} V + C_{\text{G}} \left( {V - V_{\text{G}} } \right) + C_{\text{D}} (V - V_{\text{D}} ) $$
(5)
Fig. 4
figure 4

Capacitance model for SET

where q is the charge of electron and the potential V (with Vs grounded) is.

$$ V = \frac{{C_{\text{G}} V_{\text{G}} + C_{\text{D}} V_{\text{D}} }}{{C_{\text{S}} + C_{\text{G}} + C_{\text{D}} }} + \frac{ - q\Updelta n}{{C_{\text{S}} + C_{\text{G}} + C_{\text{D}} }} $$
(6)

Now, the potential energy U of the island can be obtained by multiplying the potential by –q,

$$ U = - qV = - q\left[ {\frac{{C_{G} V_{G} + C_{D} V_{D} }}{{C_{S} + C_{G} + C_{D} }} + \frac{ - q\Updelta n}{{C_{S} + C_{G} + C_{D} }}} \right] $$
(7)

Computational procedure

In this paper, we have calculated the charging energies of IAB molecules for plotting the charge stability diagrams in both the isolated as well as SET environments, and to analyze the total energy as a function of gate voltage in SET environment. For estimating the charging energies of the molecule, total energies of the molecule for different charge states (-2, -1, 0, 1, 2) have been calculated using ATK 11.2.2 Toolkit, which based on density functional theory (DFT) and non-equilibrium Green’s function.

Here, in the computation, we have used the IAB molecule as an island in the SET, placed in between the metallic source and drain electrodes above the gate dielectric. The gate dielectric is of 3.8 Ǻ thickness and 10 ε0 dielectric constant lies on a metallic back-gate. The IAB molecule is placed 1.2 Ǻ above the gate dielectric. The separation between IAB to the drain and source electrodes is 2.8 Ǻ. The work function of electrodes is taken as 5.28 eV, which is of gold (Riviere 1966). The IAB SET architecture is shown in the Fig. 5. For the analysis of SET, zero volts source-drain bias is applied. To calculate the total energies in the SET environment, we used non-polarized DFT, and local density approximation (LDA) is taken as the exchange correlation function. Basis set double-zeta polarized (DZP) is used for wave functions expansion. Neumen boundary conditions have been applied such that the perpendicular component of the electric field is zero at the boundaries. Further, the energy zero of the molecular energy spectrum is set to absolute energy.

Fig. 5
figure 5

Impurity added benzene (IAB) based SET

Results and discussion

We have performed the analysis of benzene by adding boron and nitrogen as impurities in the gas phase (isolated) as well as in SET environment. The impurity addition has been performed in two ways, replacing the last carbon atom of benzene by boron [Fig. 6(b) (i)] and by nitrogen [Fig. 6(b) (ii)] atoms produces C5H6B (boratabenzene) and C5H6N (pyridinium), respectively. Replacing the last hydrogen atom of benzene by boron [Fig. 6(c) (i)] and by nitrogen [Fig. 6(c) (ii)] produces C6H5B and C6H5N, respectively. The charging energies for benzene and IAB molecules have been computed by calculating the total energies of the Island for different charge states of the molecule (-2, -1, 0, 1, 2) in isolated as well as SET phases and tabulated in Table 1. The charging energy is defined as energy required for charging the tunnel junction with one elementary charge. The ionization potential (EI) and electron affinity (EA) are calculated using the equations shown below

$$ E_{I} = E\left( N \right) - E\left( {N - 1} \right) = E^{0} - E^{ + 1} $$
(8)
Fig. 6
figure 6

Benzene and impurity added benzene (IAB) molecules

Table 1 Calculated charging energies of IAB molecules in isolated and SET environment

where E0 is the total energy of neutral system having N electrons, E+1 is the total energy of the positive ion having (N − 1) electrons.And,

$$ E_{A } = E\left( N \right) - E\left( {N + 1} \right) = E^{0 } - E^{ - 1} $$
(9)

Similarly, where E−1 is the total energy of the negative ion with (N + 1) electrons.

In SET phase, reduction in charging energy is due to the stabilization of the charge on the island by the electrostatic surrounding, and the reason for this stabilization is the use of dielectric material and the IAB molecule lying flat on the dielectric.

In SET, the charging energy of the island can be modified by tuning the gate voltage, which moves the energy levels of the SET in and out of the bias window. This dependence of the SET conductance on the gate voltage and the bias potential can be illustrated with help of charge stability diagram shown as Fig 7-II. These plots consists of color codes, shows the number of energy levels inside the bias window [blue (0), light blue (1), green (2), orange (3), red (4)] for a given value of the source–drain bias and the gate voltage. The conductance of SET is directly related to the number of energy levels inside the bias window. The charge stability diagrams have been plotted for all the IAB SETs. It can be seen from Fig. 7-I, that in SET environment, at zero gate voltage, C5H6N has less total energy for its positively charged state (one electron removed from C5H6N), thus more stable in its positively charged state compared to other states. C6H5B is found stable for both the neutral as well as positively charged states at zero gate potential; whereas C6H6, C5H6B, C6H5N are stable in their neutral states at zero gate potential. Figure 7(II) shows the plotted charge stability diagrams of benzene and all IAB SET. From Fig. 7(II), as expected, benzene and all the IAB SETs are not having any energy level and does not supports conduction at zero gate and source–drain bias. Due to the impurity addition, an improvement in the SET conductance has been observed in terms of gate and source–drain bias. A small change in the gate and source–drain bias does not have much importance for the conduction of benzene based SET, since the energy levels does not come into picture for small change in bias potentials for this SET; whereas, for IAB SETs (C5H6B, C5H6N, C6H5B, C6H5N based SETs), a small change in gate and source–drain potential causes the device to conduct.

Fig. 7
figure 7figure 7figure 7

(I) Total energy as a function of gate potential in SET environment. Each curve represents a charge state: blue [-2], green [-1], red [0], turquoise [1], and violet [2]. (II) Charge stability diagram in SET environment. The colors show the number of Energy levels inside the bias window for a given gate voltage. The colors represent: blue (0), light blue (1), green (2), orange (3), and red (4)

Conclusion

SET with IAB island has been modeled and compared with the pure benzene SET using the ATK 11.2.2 (Atomistix Toolkit, Quantum wise) code. In finding, we observed that at zero gate potential in SET environment, C5H6N shows more stability in its positively charged state in reference to its less total energy in its positively charged state and C6H5B shows more stability in neutral and positively charged states; whereas C6H6, C5H6B, and C6H5N show more stability in their neutral states only. The charge stability diagrams for the IAB SET have been plotted, which illustrates the dependence of conductance of SET on bias potential and gate voltage. With these findings, we conclude that the work presented here will be of great importance in understanding the role of impurity on overall performance of the SET. The IAB SETs show good improvement in the conductance in comparison to benzene based SET in terms of gate and source–drain bias.